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L INTRODUCTION 

Accurate wavelength calibration of imaging spectrometer data is essential if 
proper atmospheric transmission corrections are to be made to obtain apparent surface 
reflectance (Gao et aL, 1993). Accuracies of 0.1 nm are necessary for a 10 nm-sampling 
instrument in order to match the slopes of the deep atmospheric water vapor features that 
dominate the 0.7-23 pm region (Goetz et aL, 1991). 

The Airborne Visiblc/Infraied Imaging Spectrometer (AVIRIS) (Vane et aL, 
1993) is calibrated in the laboratory to determine the wavelength position and full-width- 
half-maximum (FWHM) response for each of the 224 channels (Chrien et al., 1990). The 
accuracies are limited by the quality of the monochromator used as a source. The 
accuracies vary from ±0.5 to ±13 nm depending on the wavelength region, in general 
decreasing with increasing wavelength. Green et al. (1990; 1992) make corrections to the 
wavelength calibrations by using the known positions of 14 atmospheric absorption 
features throughout the 0.4-23 pm wavelength region. These features, having varying 
width and intensity, were matched to the MODTRAN model (Berk et al., 1989) with a 
non-linear least squares fitting algorithm. A complete calibration was developed for all 
224 channels by interpolation. 

Instrument calibration cannot be assumed to be stable to 0.1 nm over a flight 
season given the rigors of thermal cycling and launch and landing loads. The upcoming 
sensor HYDICE (Rickard et aL, 1993) will require a means for in-flight spectral 
calibration of each scene because the calibration is both temperature and pressure 
sensitive. In addition, any sensor using a two-dimensional array has the potential for 
systematic wavelength shifts as a function of cross-track position, commonly called 
"smile" (Goetz et al., 1991). Therefore, a rapid means for calibrating complete images 
will be required. The following describes a method for determining instrument 
wavelength calibration using atmospheric absorption features that is efficient enough to 
be used for entire images on workstations. This study shows the effect of the surface 
reflectance on the calibration accuracy and the calibration history for the AVIRIS B 
spectrometer over the 1992 flight season. 

2. METHOD 

The O2 absorption feature, centered at 762 nm, was used here to test various 
methods for instrument calibration. The data used came from eight flights spanning the 
1992 flight season for AVIRIS. Since only the 762 nm band was used, only the B 
spectrometer calibration of AVIRIS was tested. 

The standard models, such as ATREM (Gao et al., 1993), used for radiometric 
correction are based on spectral atmospheric transmission values given in HT1RAN 
(Rothman et al., 1987). In ATREM the H3TRAN data are resampled to provide an equal 
wavelength spacing of 2.5 nm. Therefore, we chose to represent the O2 feature with 



2.5 nm samples, deemed adequate given the AVIRIS FWHM of 9 nm and sampling of 
approximately 10 nm (fig. 1). Three AVIRIS channels are affected by die O2 band, but 
the shape of the band is dependent on the wavelength position of each individual detector 
(fig. 2). 


The apparent shape of the O2 band is dependent not only on the individual 
detector wavelength and spectrometer resolution, but also on the shape of the spectral 
reflectance curve within the three channels affected (fig. 3). This effect has been 
recognized by others (Gao and Goetz, 1990; Carrere and Conel, 1993) and accounted for 
in the non-linear least squares fitting routines. However, non-linear least squares is a 
computing intensive technique, not suitable for processing large image data sets rapidly. 

The method chosen here makes use of the slopes of the spectral curves about the 
band;center wavelength channel (fig. 4). The model spectrum is resampled using 100 
different sets of wavelengths, each differing by 0.1 nm. Ratios of slope(left)/slope(right) 
are taken and the results for all possible ratios are shown in figure 5. The procedure for 
determining the actual wavelength of the center channel is as follows: (1) Divide the 
AVIRIS radiance spectrum by the exo-atmospheric solar irradiance spectrum. (2) 

Correct for the slope of the surface spectral reflectance by assuming a linear change in 
reflectance between the channels adjacent to the channels affected by the O2 absorption. 
(3) Calculate the ratio of the slopes of the two sides of the O2 feature. (4) Determine the 
center wavelength from the look-up table, the values of which are plotted in figure 5. ' 

3. RESULTS 

Figure 6 shows the results of an analysis of a full scene of data from the June 21 , 
1992 coverage of Blackhawk Island, Wisconsin. The histograms show the distribution of 
calculated center-channel wavelength values for data with and without surface reflectance 
corrections. The separation between the two peaks is the greatest for vegetated areas 
because the 762 nm band lies at the shoulder of the near infrared (NIR) plateau in 
vegetation reflectance. The width of the distributions can be attributed to system noise 
and variations in surface materials. Therefore, it is important to calculate the band center 
for a large number of pixels to get an accurate wavelength value. Figure 7 shows the 
results calculated from vegetated and unvegetated areas within the Blackhawk Island 
image. The offset between the two curves is 0.4 nm, and can be attributed to the 
assumption of a linear surface reflectance curve over the 50 nm wavelength region 
bridging the O2 absorption feature. The width of each distribution can be attributed to 
system noise. 

A number of images from the 1992 AVIRIS flight season were processed. The 
center wavelength position was constant for scenes from two different overflights taken 
on the same day. Figure 8 shows the variation of the center wavelength position of 
channel 42 with time. The wavelength calibration file accompanying the image data 
gives a value of 764.0 nm. The maximum variation is only 3 nm, but this becomes 
significant when atmospheric correction is attempted. Figure 9 shows the results of the 
application of the ATREM program to data from Gainesville, Florida taken July 8 under 
very high humidity conditions. The proper calibration is most important here because of 
the great depth of the unsaturated atmospheric water absorption features. The properly 
calibrated data set has far fewer artifacts than the one using the published wavelength 
values. AVIRIS contains an onboard calibration capability using filters. The data have 
not yet been checked to confirm the shift. 

4. CONCLUSIONS 

We have shown that, at least for the 762 nm O2 absorption band, it is possible to 
calibrate AVIRIS , using in-flight data, to an accuracy of ±0.2 nm regardless of the 
surface reflectance. This technique could be extended to other atmospheric features in 
the wavelength range covered by AVIRIS and HYDICE. The rapid analysis made 



possible by our algorithm will make it feasible to analyze each image from sensors such 
as HYDICE, for which systematic wavelength calibration changes during flight are 
anticipated. 
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Figure 1. The oxygen feature at 
H1TRAN and AVIRIS resolution and 
sampling. 
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Figure 2. Oxygen band with sensor 
center wavelengths varying by 2 nm. 
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Figure 3. Oxygen transmittance and 
vegetation reflectance spectra. Note 
the increasing slope in the vegetation 
spectrum. 
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Figure 4. The left and right sides of 
the oxygen feature are used in 
calculating the ratio. 
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Figure 5. Ratios of two sides of the 
oxygen feature for 100 different 
sensor wavelength calibrations. 
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Figure 6. Oxygen band center 
histograms for the full Blackhawk 
Island scene, corrected and uncorrected 
for surface reflectance. 
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Figure 7. Histograms of the oxygen 
band centers for vegetated and 
nonvegetated areas near Blackhawk 
Island. 
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Figure 8. Calculated oxygen band 
center wavelengths for eight scenes 
over the 1992 flight season. 



Wavelength (nm) 


Figure 9. Reflectance spectra in 
AVIRIS spectrometer B calculated 
using ATREM from Gainesville, 
Florida. Wavelengths for the corrected 
spectrum were shifted by 2.8 nm. 
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ABSTRACT 

A method is described for rapid analysis and correction of AVIRIS and HYDICE data for 
spectral shifts induced by the flight environment on the sensor. The method makes use of 
the O 2 absorption at 762 nm but can be extended to other features elsewhere in the 0.4- 
2.5 pm region. The method yields a precision of ±0.2 nm. Data from the AVIRIS 1992 
flight season were analyzed and show drifts of up to 3 nm in the B spectrometer. 


